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ABSTRACT 

We derive the fraction of blue galaxies in a sample of clusters at z < 0.11 and 
the general field at the same redshift. The value of the blue fraction is observed to 
depend on the luminosity limit adopted, cluster-centric radius and, more generally, 
local galaxy density, but it does not depend on cluster properties. Changes in the blue 
fraction are due to variations in the relative proportions of red and blue galaxies but 
the star formation rate for these two galaxy groups remains unchanged. Our results 
are most consistent with a model where the star formation rate declines rapidly and 
the blue galaxies tend to be dwarfs and do not favour mechanisms where the Butcher- 
Oemler effect is caused by processes specific to the cluster environment. 
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1 INTRODUCTION 
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The observation bv lButcher fc Oemlerl lll984) that clusters 
of galaxies contain a larger fraction of blue galaxies at pro- 
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gressively higher redshift (the Butcher-Oemler effect) con- 
tributed to the establishment of the current view of clusters 
as sites of active galaxy evolution, in which environmental 
influences dramatically alter the morphologies and star for- 
mation histories of their members. 

However, as the number of clusters observed has in- 
creased, it has become clear that the Butcher-Oemler effect 
is not solely or simply an evolutionary trend. The large scat- 
ter observed in the blue fractions for clusters in narrow red- 
shift ranges l| Smail et alJ ll993 : iMargoniner &c de Carvalhd 
I2OOOI : iGoto et^dTl2003air implies the existence of environ- 
mental effects, which may compete with, and possibly 
mimic, evolutionary trends, if the cluster samples evolve 
with redshift. 

Studies of the ori ginal Butcher-Oemler sa mple lend sup- 
port to this scenario: iNewberrv et all (^8^ measured the 
surface densities and velocity dispersions of seven Butcher- 
Oemler cluste rs and found that t he hig h and low redshift 
sample differ. lAndreon fc Ettoril il999ll show that the X- 
ray luminosity of high redshift Butcher-Oemler clusters is 
higher than that of low redshift clusters, although the X-ray 
luminosity function of clusters does not evolve to z ~ 0.8 
jHolden et al.ll2002h . arguing that the high redshift sample 
does not repre sent the progenitors of modern-day clusters . 
Similarly, both ISmail ~ t al. (199^ and lFairlev et alJ ll2002ll 
find a low blue fraction in their samples of X-ray selected 
clusters. 

The blue fraction has also been observed to depend 
on a number of other factors, such as: the luminos ity limit 
used and the clust er centric distance (Ellingson ct alj|200ll 
iGoto et al1l2003af ). richness jMarg onincr ^ al.ll2o6lf) . clus- 
ter concentration iButcher fc OemleniigslT and, p ossibly, 
the presence of substructure jMetevie^^l] lioOOf) . These 
findings point to the necessity of understanding environ- 
mental effects on the blue fraction in order to disentangle 
evolutionary trends from selection biases in studies of clus- 
ters at high redshift. 

The aim of the present paper is to study how the frac- 
tion of blue galaxies varies as a function of cluster properties 
and in the field at the same redshift in the local universe, 
in order to analyse the effects of environment in isolation 
from the evol utionary trends tha t give rise to the B u tcher - 
Oemler effect. iLewis et all J2002l) and iBaloeh et alJ (l2003l) 
have recently discussed how the star formation rate, as mea- 
sured from the Ha equivalent width, varies with environ- 
ment within the 2dFGRS sample. However, colours provide 
a measure of the average star formation histories over longer 
timescales than Ha and may therefore better reflect the in- 
fluence of environment on galaxy properties, especially if 
the mechanisms responsible for the Butcher-Oemler effect 
operate on longer t ime frames, as is the case for harassment 
jMoore et alJll99et) . for instance. 

The structure of this paper is as follows: in the next 
section we present our analysis; selection of clusters, clus- 
ter members and a comparison field sample, definition of 
the radial and luminosity limits used, fc- correct ions, colour- 
magnitude relations, density measurements and a discus- 
sion of how the blue fraction and its error are measured. 
The main results are detailed in §3. below and we dis- 
cuss these in §4. Throughout this paper we adopt a cos- 
mology with Qm = 0.3, f^A = 0.7. It is normal to define 
h = //o /100 km s~^ Mpc~^: here we suppress the h scaling. 



so that the full explicit meaning of Mpc in length units is h 
Mpc; for absolute magnitudes M stands for M + 5 logj^Q h. 



2 ANALYSIS 

In this work, we analyse how the blue fraction of galaxies 
varies as a function of cluster properties and field density 
at z < 0.11; our purpose is to provide a local reference for 
studies of galaxy colour evolution in high redshift clusters 
(the Butcher-Oemler effect). It is therefore useful to repro- 
duce, in our analysis, so me of the features of the original 
iButcher fc Oemlen (1198411 study, as this is the basis for most 
current work on the evolution of stellar populations of clus- 
ter galaxies. 

The original analysis of iButcher fc Oemleil il984|) de- 
fined blue galaxies as being (i) within a radius contain- 
ing 30% of the cluster population; (ii) brighter than a no- 
evolution fc— corrected Mv = —20 and (iii) bluer by 0.2 mag- 
nitudes in B — V (no-evolution and fc— corrected) than the 
colour-magnitude relation defined by the cluster early-type 
galaxies. Below we describe how we implemented these pre- 
scriptions for our dataset. 

2.1 Cluster and member selection 

The clusters a n alyzed here are those studied by 
iDe Propris et al.l tOQSi) in their paper on the com- 
posite galaxy luminosity function. This sample consists of 
60 clusters at 2; < 0.11 containing at least 40 spectroscopic 
members and with average 85% completeness. Although 
these clusters share, to some extent, the biases of the 
Abell, APM and EDCC catalogues from which they were 
originally drawn (see De Propris et al. 2002 for details), 
these 60 objects provide an at least approximately complete 
and volume limited ensemble of nearby clusters, spanning a 
large range of properties (such as richness and velocity dis- 
persi on) determined from the same data (IDe Propris et alJ 

Cluster membership was determined via a 'dou- 
ble gapping' method, as described in greater detail in 
iDe Propris et all i2002h : for each putative cluster centre and 
redshift, we first required that the likely cluster members be 
surrounded by 1000 km s"'^ gaps on either side of the red- 
shift distribution; we next computed a first estimate for the 
mean cz and velocity dispersion (ar) and ranked all galax- 
ies in order of distance from the mean cz. If galaxies had a 
'gap' in redshift space from their neighbour larger than our 
first estimate of ar all galaxies with a faster or slower veloc- 
ity were excluded. This procedure yields a sample of cluster 
members which is likely to be relatively uncontaminated by 
interlopers (i.e. field galaxies with appropriate redshift but 
not dynamically bound to the cluster). 

2.2 Aperture and luminosity selection 

We calculate rao (the original aperture used by Butcher & 
Oemler 1984) by using all cluster members within a 3 Mpc 
radius from the central galaxy. We then choose an aper- 
ture than includes 30% of these galaxies. It is observed that 
the value of rao varies depending on cluster dynamics and 
central concentration. For this reason it may be generally 
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more appropriate to use apertures based on cluster dynami- 
cal properties, in order to insure that all gal axies considered 
share the same environment in all clusters. lEUingson et al.l 
l|200ll) advocate using r2oo, the radius at w hich the cluster 
is 20 times denser than the general field iCarlberg et alJ 
119971) . or its multiples. This is defined as: 



''200 



WH{z) 



(1) 



where H{z) represents the redshift-dependent Hubble con- 
stant. In the following, we will use both r2oo and r2oo/2 for 

our analysis. 

The original limit by iButcher fc Oemleil il984l) cor- 
responds to approximately 1.8 magnitudes below the M* 
point. Here, we use a magnitude limit 1.5 magnitudes 
fainter relative to our computed M* and we also analyse 
the effects of using fainter limits. This allows us to study a 
more consistent sample of galaxies in all clusters. 



2.3 fc-corrections 

Bj — Rf colours for all 2dFGRS galax ies were derived usin g 
data from the SuperCosmos survey iHamblv et al.ll200ll) . 
Given the colours and redshifts we derived a fc-correction in 
the following way. 

The KCORRECT package iBlanton et alj l20'03a^ fits a 
combination of templates to galaxy colour data, and thus 
yields a consistent fc-correction. Although only two bands 
are available for the 2dFGRS, this approach was validated 
by fitting the full DRl ugriz data, then comparing with the 
result of fitting gr only. The differences in K-correction are 
generally at the 0.01 mag. level. 

The main problem with KCORRECT is that the tem- 
plates are unable to describe extremely red galaxies {Bj — 
Rf = 1.26 corresponding to g — r — 0.87), and so a 
consistent fc-correction cannot be obtained. In such cases, 
we took an alternative a pproach, based on the models of 
iBruzual fc Charlod ll2003l) . Each galaxy was modelled by a 
single burst, varying the age until the observed colour was 
matched at the given redshift. This age is degenerate with 
metallicity; in practice, we assumed 0.4 times Solar metallic- 
ity {Z = 0.008), except for very red galaxies where a current 
age of > 13 Gyr was required. In these cases, the current 
age was set at 13 Gyr and the metallicity increased until the 
colour was matched. For red galaxies, fc-corrections deduced 
in this way agree almost exactly with the results of KCOR- 
RECT, and we were able to match smoothly from one to the 
other to deal with the very red galaxies that KCORRECT can- 
not model. The results differ for extremely blue galaxies (at 
the 0.1 mag. level), but a single burst is a rather unrealistic 
model in such cases and the results of KCORRECT are to be 
preferred where they are consistent. 

The fitting functions given below describe the k- 
correction results for all colours of practical interest at 
z < 0.3, to within a maximum error of about 0.02 mag. 
Given the intrinsic uncertainties in both photometric mod- 
elling and calibration, and in the interests of clarity, these 
residuals were ignored and the fit was treated as exact. 



k{Bj) = (-1.63 + 4.53x}y ~ (4.03 + 2.01a;)y 

z 



where 



l+(10z)4 

k{RF) = (-0.08 + lA5x)y - (2.88 + 0.48x)y'^ 

X = Bj ~ Rp 
y = z/{l + z) 



(2) 



(3) 



(4) 



(5) 



2.4 Colour-magnitude relations 

For all clusters, we determine the colour-magintude relation 
by carrying out a least absolute deviation regression fit t o 
the observed colour distribution jArmstrong &: KunJl97^ . 
The average observed slope is about —0.019 ± 0.003, con- 
sistent with the relation observed in the Coma galaxies for 
V — R (Eisenhardt et al. 2004, in preparation). There ap- 
pears to be some cluster-to-cluster variation in the slope, but 
this may be due to small number statistics in some cases. 

We carry out the same analysis for field galaxies, but we 
divide the sample according to local density (see below for 
details). The slope for the field galaxies is ~ —0.013± 0.001. 
These slopes were calculated robustly and jackknife resam- 
pling was then used to compute the errors on the values of 
the slopes. The slopes of the relations for field and cluster 
galaxies appear to be different at about the 2a level, but we 
do not regard this as compelling (cf., Hogg et al. 2004 for a 
similar study from Sloan colours). 

The 2dFGRS SuperCosmos colours can be transformed 
to the SDSS g and r bands using the following equations: 



Bj = ff-h 0.15 + 0.13(5 -r) 



Rf 



0.13 



(6) 



(7) 



A difference in i? — 1/ of 0.2 magnitudes from the ridge line 
defined by early- type g alaxies is equivalen t to a 0.2 magni- 
tude difference in g — r llGoto et alJl2003a^ . From the above 
equation we then derive that the corresponding difference 
in Bj — _Rf is 0.23 magnitudes. Fig. 1 shows the colour- 
magnitude relation and the colour distribution, marginalised 
over the derived colour magnitude relation, for a few repre- 
sentative clusters. 



2.5 Calculation of local density 

For field galaxies, the only relevant 'environmental' prop- 
erty is their local density. We calculate this by computing 
the number of galaxies to Mb, = —19 in an 8 Mpc sphere 
centred on each galaxy, with appropriate completeness cor- 
re ctions. A detailed d escription of this procedure is given 
in lCroton et alJ (|20o3). Unfortunately, this assumes that all 
radial velocities are due to the smooth Hubble flow and can 
be used as proxies for distance. This is not applicable in the 
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Figure 1. Colour-magnitude relation (left) and histogram (right) of marginalised colour distribution (opposite each colour-magnitude 
relation) for four representative clusters. The arrow in the histogram indicates the blue/red separation (see text) 



cluster environment, where the galaxies have turned around 
and no longer participate in the cosmological expansion (see, 
e.g., Cooperstock et al. 1998). For clusters, therefore, we as- 
sume that all members are effectively within r2oo and calcu- 
late the density in an 8 Mpc sphere containing only cluster 
members. This places the density for cluster galaxies on the 
same scale as that used for field galaxies. As we note below, 
this is likely to be a slight underestimate of the actual den- 
sity (as some members may exist beyond the virial radius). 



2.6 Calculation of the blue fraction 

Calculation of the error in the derived blue fraction has been 
carried out with a variety of recipes. Here we derive the ap- 
propriate formulation for our sample of cluster members. 
Note that if the numbers of red and blue galaxies are deter- 
mined via background subtraction, extra contributions to 
the error statistics due to clustering need to be included. 

The blue fraction is defined as the ratio of m blue galax- 
ies observed out of n total galaxies. Assuming that m and 
n obey Poisson statistics, the blue fraction is: 



(8) 



and its likelihood has the same form independent of whether 
we assume Poisson or binomial statistics (i.e. with n fixed 
in advance): 



(9) 



whose maximum is, trivially, m/n. Therefore the variance 
of the blue fraction is: 



(f\nL 



m{n — m) 



(10) 



This is incorrect for m — because dL/df 7^ for m = 0. 
In that case we ask what yields likelihood equal to 
exp(— 1/2 fmax)- This is approximately l/2n, which is a rea- 
sonable error bar to adopt for the m = case. 

In our analysis we only use spectroscopically confirmed 
cluster members. Since the redshift sample is Bj selected, 
without regard to colour, this should not bias the result, 
especially given our high completeness (cf. Elhngson et al. 
2001 for a similar approach to data in the CNOC2 survey). 



3 RESULTS 

3.1 Luminosity and radial dependence 

The blue fraction depends on the luminosity limit 
and on the aperture used to include member s 
jMargoniner fc de Carvalhl 1200(1 : lEUineson et ai\ I2OOII ). 
We analyse these dependencies here, using our sample. 
Fig. 2 shows the blue fraction, calculated to M* + 1.5, vs. 
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Figure 2. The blue fraction to A/* + 1.5 for all clusters, in all 
three apertures considered 




Figure 3. Mean values of the blue fraction as a function of 
luminosity and aperture 



z for all three radii we consider: rso, 1200 and 1200 /2 (where 
M* is taken from De Propris et aL 2003a). This plot shows 
that there is no dependence of the blue fraction on redshift 
and that therefore our sample of objects is adequate to 
study environmental trends in isolation from the effects of 
evolution. Unless otherwise specified, we refer to this blue 
fraction as jh for the remainder of this paper. 

Inspection of Fig. 2 shows that there is consider- 
able scatter in the derived blue fraction. This is consis- 
tent with previous studies of local and X-ray selected sam- 
ples llSmail et al.lll998l: iMargoniner et al.ll200l] : iGoto et alJ 
l2003a^ . We compare the distribution of points with a Gaus- 
sian distribution with mean and standard deviation derived 
from a weighted fit to the data. The distribution shows sig- 
nificant skewness and a ^ test yields probabilities below 
10^'^ that the data points are consistent with a single value 
of /;,. This suggests that the scatter observed is a real prop- 
erty of the sample. 

We plot the mean values of fj, as a function of the 
magnitude limit used for all three apertures in Fig. 3. A 
t— test shows that the distributions have significantly differ- 
ent means: the blue fraction increases as a function of radius 
and as a functi on of luminosity. This is c onsi stent with previ- 
ous fin dings bv lMargoniner et al.l (|20^) and lEUingson et alJ 
l|200lfl and implies that the the blue galaxies are intrinsi- 
cally faint objects and reside preferentially in the cluster 
outskirts. 



3.2 Dependence on cluster properties 

We will next consider how the blue fraction depends 
on cluster properties derived from our original study 
JPe Propris et all2002l) : Bautz-Morgan type, richness, pres- 
ence of substructure, velocity dispersion and concentration. 

The Bautz-Morgan type measures the dominance of the 
brightest cluster galaxy relative to the rest of the cluster 
population. It is used as an indicator of dynamical evolu- 
tion, where the brightest galaxy evolves via the merger of 
M* galaxies at the cluster center. Velocity dispersion pro- 
vides a measure of the cluster mass (assuming the clusters 
to be virialized) and of the relative speed of galaxy encoun- 
ters (which is important for some of the mechanisms for 
the origin of the Butcher-Oemler effect). Richness, as mea- 
sured from the number of galaxies brighter than Mi,j = — 19 
iPe Propris et alJl2003ah . is an indicator of cluster mass 
and mean density. We consider the probability that a clus- 
ter contains substructure, using the Lee-Fitchett statistics 
(Fitchett 1988). Although there is no metric to determine 
the amount of substructure a cluster contains, the probabil- 
ity of its being significantly substructured provides a mea- 
sure of the amount of recent merging a cluster has under- 
gone and therefore may allow us to estimate how recently 
the cluster has formed. Finally, we consider the dependence 
on cluster concentration, defined as 

C = Iog(r6o/r2o) (H) 
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Figure 4. Dependence of the measured blue fraction on cluster 
properties for tlie r200 /2 aperture 



where rgo and reo are the radii containing 20% and 60% of 
the cluster population, respectively. 

Fig. 4 shows how the blue fraction (to M* + 1.5) depends 
on these five variables. We only plot the case corresponding 
to r2oo /2 for economy of space, but the results are similar in 
all three apertures. In no case is there any clear evidence of 
any dependence of the blue fraction on Bautz-Morgan type, 
richness, velocity dispersion, probability of substructure or 
concentration. This is confirmed by a non-parametric Spear- 
man test, which returns low correlation probabilities in all 
instances. 



3.3 Density dependence in the field 

We have selected a sample of field galaxies from the 2dFGRS 
galaxy redshift survey with the same redshift distribution as 
the clusters. For each galaxy we derived the local density in 
an 8 Mpc sphere (Croton et al. 2004). From this sample we 
determine the blue fraction as a function of local density in 
Fig. 5. 

The blue fraction follows a well defined linear trend, 
becoming higher in low density regions. This trend can be 
extended to clusters, whose behaviour is then simply the 
more extreme version of the relation between ft and local 
density as observed in the general field. 

Galaxies then display a bimodal distribution in colour 
space, with well-defined red and blue wodges, whose relative 
populations change with density. This represents a confirma- 
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log p (galaxies Mpc " ) 

Figure 5. The blue fraction in the general field as a function of 
local density (filled circles) and the mean for all clusters within 
r2oo/2 (open square). Density is defined within 8 Mpc spheres as 
in Croton et al. (2004). For clusters we use the number of galaxies 
within r200 and assume that all cluster members are within this 
sphere in our calculation of density. Because this is likely to be a 
slight underestimate, we indicate in the figure that this density is 
only a lower limit. 



tion of the trends observed by the SPSS in their more ext en- 
sive photometry jHogg et al.ll2003l: iBlanton et al.ll2003bl) . 

It is interesting to compare how the relative star for- 
mation rates, as measured from Ha equivalent widths, vary 
with density or with the only properties that, in clusters, 
appear to affect the blue fraction, luminosity and radius 
(computing absolute star formation rates is difficult with 
2dF data, as the spectra are not easy to flux calibrate, but 
Ha equivalent widths can be used to provide an estimate of 
the relative star formation rate within the sample). The Ha 
equivalent widths were computed using Gaussian line fitting 
with a small (2 A) correction for underlying line absorption. 
No dust extinction correction was carried out . A full de- 
scription of the procedures used can be found in lLewis et alJ 
(2002), especially their section 2.4. 

The variation of Ha equivalent width with density, lu- 
minosity and radius is shown in Fig. 6 and 7. for red and 
blue galaxies in both environments. We see that, while the 
blue fraction changes the mean relative star formation rate 
does not. The bimodal behaviour observed in t he colours ex- 
tends to the Ha equivalent width, as noticed bv lBaloeh et alJ 
i2003h . The changes in blue fraction and Ha equivalent 
width as a function of density are due to changes in the 
relative fractions of quiescent and star forming galaxies. 



4 DISCUSSION 

We have determined the fraction of blue galaxies in a sam- 
ple of nearby clusters and in the general field at the same 
redshift. Although the blue fraction varies considerably from 
cluster to cluster, we find that the variation is a real property 
of the sample and not due to statistical noise. 

The mean value of the blue fraction is 0.11 ± 0.01 
for galaxies within rso and 0.13 ± 0.01 for galaxies 
within r2o o/2. These are higher than the v alue of 0.03 ± 
0.09 from iMargoniner fc de Carvalhol J2000[). ~ 0.07 fi-pm 
iMargoniner et alJ J200lll and ~ 0.07 from Ipimbblet et alJ 
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Figure 6. Median Ha equivalent widths (in Angstroms) for blue 
and red galaxies in the field as a function of density. Errors are 
interquartile means. The errors for the red galaxies are smaller 
than the symbols. 



I|2002fl at 2 < 0.11, although they are within the errors. One 
possibility for this discrepancy lies in the possibility of blue 
interlopers from the field population being included in our 
redshift sample jPiaferio et ahlliooj) . We have calculated 
the level of contamination by integrating the field luminos- 
ity function of Madgwi ck ct al. (2002i) over the appropriate 
luminosity range and within the pseudo-volume defined by 
the aperture used and the range in velocities spanned by 
cluster galaxies. The blue fraction for this sample was cal- 
culated from our density-dependent blue fraction. The av- 
erage level of contamination is 0.08 ± 0.06 but this depends 
critically on assumptions concerning the mean density of 
galaxies in cluster outskirts (which influence both the nor- 
malization and the fraction of blue interlopers as per Fig. 5). 
In addition, there are contributions to the error from both 
Poisson and clustering statistics. The number of blue galax- 
ies is, in any case, small and therefore the level of contami- 
nation is uncertain and the error in its estimate quite large. 
We therefore ignore the issue of possible contamination of 
our spectroscopic sample of cluster members. 

• A first conclusion to be drawn from the derived blue 
fractions is that few clusters have > 40% and none has 
fb > 60% . 

Since this sample of clusters is, at least to first order, 
complete and volume limited, it represents a fair sampling 
of the range of blue fractions encountered in cluster envi- 
ronments. By contrast, samples of clusters at high redshift 
contain at least a few objects with blue fractions in excess 
of 40% (Fairley et al. 2002, La Barbera et al. 2003 and 
references therein). This would suggest that the observed 
evolution is real: contamination in our sample (if real) 
should only make the local blue fraction lower, while high 
redshift clusters are drawn from the high richness envelope 
that contains few local clusters with high /;,. However, it 
is possible that optically selected clusters at high redshift 
tend to contain higher fractions of blue galaxies because 
these make them more conspicuous in blue plates. 

• The blue fraction appears to depend on the lumi- 
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Figure 7. Median Ha equivalent widths (as in Fig. 6) for blue 
and red galaxies in clusters as a function of luminosity and radius. 
Errors are interquartile means. 



nosity limit and cl uster centric radius used ( Fig 



^ 3)> as 

previously found bv [ Margoniner fc de Carvalhd ll2000fl and 
lEUineson et al.1 J200ll) ^ 

This behaviour is similar to that observe d for dwarf 
galaxies in Coma and Abell 2218: iPracv et al.l |2003) have 
shown that dwarfs have a steeper luminosity function at 
larger cluster centric radii and are preferentially found at 
large distances from the cluster core, with these trends 
being stronger for the faintest dwarfs. Similar trends 
may have been observed f or galaxies in the Coma cluster 
jBeiiersbergen et alJ l2002ll and appear to persist in the 
z ~ 0.4 CNOC sample. To the extent that present-day 
clusters and their populations represent proxies for high 
redshift objects, these observations imply that a proportion 
of the blue galaxies are intrinsically low luminosity objects. 



• There is no dependence of the blue fraction on clus- 
ter properties (Fig. 4). The lack of correlation of the blue 
fraction with cluster properties suggests that the star for- 
mation rate, or its decrease, is not related to the large scale 
structure: therefore, mechanisms which involve cluster-wide 
effects, such as ram stripping by the cluster gas, tides, which 
depend on the cluster mass, or harassment, whose efficiency 
is proportional to velocity dispersion, are disfavoured by the 
present study, as well as other explanations that depend on 
processes specific t o the cluster environment. However, as 
iBalogh et alJ (|20o3) show, the star formation rate is affected 
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by local processes and close interactions are a viable mech- 
anism for the origin of the blue fraction. 

The above is ap parently in contrast with some previous 
studies: iMetevieret al. (2000,) claim that clusters containing 
large amounts of substructure have higher blue fractions 
and interpret this in terms of a shock model during colli- 
sions of subclusters. However, their result is based on three 
clusters and may not be vali d for the general po pulation. 
iMargoniner et al.l teOOlfl and iGoto et al.l ijiooHi) suggest 
that the blue fraction depends on richness. One possibility 
is that, since they use a fixed 0.7 Mpc aperture, and given 
the radial dependence observed in Fig. 3, that a spurious 
richness dependence may be induced by aperture effects. 
We have derived blue fractions for a 0.7 Mpc aperture in 
our clusters and tested for correlation with richness but 
failed to detect a significant signal. Similarly, we failed to 
detect any correlat ion with richness if we use the same 
luminosity range as IMargoniner et al.l (1200111 . However, we 
see in Fig. 4 that rich clusters tend not to have large blue 
fractions and we also observe that the blue fraction in the 
0.7 M pc aperture shows large scatter. If the sample used by 
IMargo niner et al. ( JOoJ) and iGoto et al. ( 20P3a) is biased 
towards richer clusters at high redshift, the combination of 
these two effects may produce a spurious correlation. 

• We observe that the blue fraction exhibits a strong de- 
pendence with local density in the general field (Fig. 5) and 
that the cluster value represents a continuation of this trend 
to higher density regimes. 

This suggests that the same processes are responsible 
for the Butcher-Oemler effect in all environments and that 
these mechanisms vary smoothly as a function of density. 
Again, this implies that cluster-specific processes are not 
likely causes of changes in the blue fraction, but local 
mechanisms whose effi ciency varies smoothly with density, 
such as interactions jLavery fc Henrvl Il988l : ICouch et alJ 
Il998l) . may be viable explanations. 

• The relative star formation rate for the blue galaxies 
does not appear to decrease as a function of density (Fig. 6) 
or as a function of luminosity or radius in clusters (Fig. 7) . 

This implies that changes in the blue fraction are solely 
due to changes in the relative fractions of red (quiescent) 
and blue (star-forming) galaxies. The trends we observe 
can be roughly reprodu ced by using the den sity-dependent 
luminosity functions of ICroton et al.l (12004) for early- and 
late-type galaxies and the typ e-dependent cluster lumi- 
nosity functions of iDe Propris e t al. (2003a) and simply 
assuming that early-type galaxies are red and late types 
ar e blue. This suggests th at the simple model presented 
in iDe Propris et al.l (l2003a^ . where galaxies simply moved 
between spectral types without number or luminosity 
evolution, provides, heuristically, a good representation of 
galaxy evolution. 

One possible caveat is that, by selecting blue galaxies, 
we have automatically selected for galaxies with Ha emis- 
sion, while morphologically selected samples s how declining 
star formation rates as a function of density dGomez et alJ 
l2003l:lGoto et al However, the observation that blue 

cluster galaxies have the same Ha equivalent width as their 
counterparts in the field is non-trivial, because lower (but 



non-zero) star formation rates will lead to weaker Ha emis- 
sion, while the colour would remain blue by our definition. 
This is in contrast, for instance, with models where galaxies 
are 'choked' in clusters (Balogh et al. 2000). 

The above argues for a model where star formation 
declines over relatively short timescal es, leading to the b i- 
modal distribution in Ha observed bv lBalogh et al.l i2003l) . 
and galaxy colours evolve quickly to the red envelope, 
producing the colo ur bimodality observed here and in 
the SPSS (iHogg et al. 2003 : Blanton et al. 2o63bl) . As in 
iDe Propris et al.l (l2003al) this is possible if the optical colour 
is dominated by the young population but the majority of 
light is provided by the underlying old stars, so that once 
the star formation is extinguished ga laxy co lours quickly 
evolve on to the passive locus. Sh iova et alJ (^002) have 
shown that truncating the star formation of spiral galaxies 
in the field leads them on to the colour-magnitude relation 
defined by the Coma E/SO galaxies, but that this process is 
inefficient for the more massive galaxies and is viable only 
for low mass spirals. When we consider the radial and lu- 
minosity trends, one possible interpretation is that a large 
fraction of the blue galaxies are actually dwarfs undergoing 
episodes of star format ion, as also suggested by a number of 
other lines of evid ence (Rako s et alll99'^ICouch et alJl998l : 
iDe Propris et al']|2bo3b) . However, we caution the reader 
that since our sample of cluster and field galaxies is local, 
we cannot properly discuss the origin of the blue fraction at 
higher redshift. 

A complication in interpreting these data is that the 
observed correlations of and Ha with density on large 
scales imply that the environment today is not affecting the 
properties of the population and therefore we are unable 
to observe galaxy evolution in progress but just its end re- 
sult. This implies that identifying the mechanisms respon- 
sible for the Butcher-Oemler effect in the local universe is 
problematic. Surveys at intermediate redshifts (e.g. DEEP2, 
VIMOS) may be able to witness the main phases of galaxy 
evolution. 
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